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Abstract. The dynamics of atom-field entanglement for a system composed of two 
atoms resonantly coupled to a single mode leaky cavity field has been investigated 
beyond rotating wave approximation (RWA). By using monogamic relation for 
entanglement of formation (EOF) as well as the lower bound of EOF for bipartite 
mixed states in higher dimensions, contrary to the RWA case, the atom-field system 
in the steady states is found to be entangled in the strong-coupling regime and the 
entanglement can grow as a function of atom-field coupling constant. 



1. Introduction 

Rabi model describes the interaction of a two-level atom (a qubit) with a quantized 
single mode cavity field (a harmonic oscillator) [1]. It has been recently used to 
describe many physical phenomena, such as an electron coupled to a phonon mode [2] , a 
superconducting qubit interacting with a nanomechanical resonator [3J, a transmission 
line resonator [1] or an LC circuit [5], etc [6]. Due to the excitation non-conserving 
nature of the Rabi model, it is hard to obtain an elegant formula for the eigen-spectrum 
of the Rabi model; only recently Braak was given the spectrum in terms of a solution of 
a transcendental equation with a single variable [7]. For the weak atom-field coupling 
regime and the nearly resonant condition, it is legitimate to drop the excitation non- 
conserving terms, the so-called counter-rotating terms, in the atom-field interaction. 
This approximation is known as rotating wave approximation (RWA) and has been 
well applied recently in most optical settings both experimentally and theoretically. 
Under RWA, Rabi model reduces to Jaynes-Cummings model [8] which is exactly 
diagonalizable. On the other hand, with the recent technological advances in the circuit 
and cavity QED systems [21 HI [5|, [6] , it is now possible to engineer the strong atom- field 
couplings which makes crucial to consider the effects of counter-rotating terms on the 
dynamics of atom-field systems. 

In the present study, we consider a model of two noninteracting qubits strongly 
coupled to a harmonic oscillator in a dissipative cavity. Our main focus is to investigate 
the effects of counter-rotating terms on the dissipative dynamics of entanglement 
between an atom and the cavity mode by using the recently derived Lindblad type 
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quantum optical master equation P] which is adapted to the strong-couphng regime 
under Born and Markov approximations. The entanglement between atoms in the same 
setting has been recently investigated by many groups in the strong-coupling regime 
by using different master equation approaches [TOl [TTl [121 [131 [IH [HI [H]. It has been 
shown that the entanglement between the atoms displays richer dynamical features in 
comparison to the RWA one. In particular, counter-rotating terms are found to lead 
to steady state entanglement between the atoms even for the initial states (including 
zero-excitation state) that have no overlap with the decoherence free subspace of the 
system [10]. Different from the previous studies [101 [HI [I2l [111 [IH [151 [16], we focus on 
the entanglement between atom-field system in the strong-coupling regime. Under RWA 
conditions, even the entanglement between the atoms can reach a high non-zero value in 
the long-time limit for certain initial states, the cavity mode always becomes separable 
from the atoms in the steady states. By using monogamic relation for entanglement of 
formation (EOF) [T71 [18] as well as the lower bound of EOF for bipartite mixed states 
in higher dimensions [I9] , we have shown that contrary to the RWA case the atom- field 
system in the steady states can become entangled in the strong-coupling regime and the 
entanglement can grow as a function of atom- field coupling constant. 

The paper is structured as follows. In Sec [21 we introduce the model. The results 
are presented in Sec. [31 We conclude with the main finding of the paper in Sec. [H 

2. The Model 

The Rabi Hamiltonian that describes the strong interaction of two identical qubits (A 
and B) with a quantized cavity field F can be written as (/I = 1) [1]: 

HR=^Ua ^ ai+Ufa^a + g ^ ai{a + a'^), (1) 

j=A,B j=A,B 

where Ua (ujf) is the qubit transition (cavity field) frequency, a (a''") is the 
annihilation (creation) operator of the cavity field, (o"^) is the Pauli matrix in the 
z (x) direction and g is the coupling constant between the atom and the cavity field. The 
terms, aia, a]_a\ in the interaction part of the Hamiltonian are called counter- rotating 
terms (CRTs) and they do not conserve the excitation number; even the atoms and the 
field are initially in their ground states, the time evolution of the mean excitation number 
is nonzero [H]. On the other hand, in the weak-coupling regime where g << Uf, Ua, the 
effects of CRTs are negligible and these terms can be dropped form the Hamiltonian. 
This approximation is known as rotating wave approximation (RWA) which leads to 
two-qubit Jaynes-Cummings Hamiltonian [8j: 

Hjc = ^^^a ^ ai + Ufa^a + g ^ [a{a + ala^] . (2) 

j=A,B j=A,B 

The photon loss due to the imperfections of the cavity mirrors should be taken into 
account and it is generally represented by the " standard^ Lindblad master equation: 

P = -ilHjc, p] + K,D[a]p, (3) 
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where D[m\p = | (2mpm^ — pim)m — m^mp^, p is the atom-field density matrix and k 
is the photon leakage rate. In the weak-coupling regime where RWA can be applied, this 
master equation is well studied and can be used accurately to describe many cavity and 
circuit QED systems. Contrary to the numerous investigations of this master equation in 
the strong-coupling regime in an ad hoc manner (see Ref. [10] and references therein), it 
was recently shown that it is at the root of unphysical effects [9]; it creates excess photons 
in atom-field system and produces spurious qubit fiipping. The main reason of these 
unphysical predictions obtained from standard Lindblad master equation in the strong- 
coupling regime is to neglect the qubit-field interaction while deriving the dissipator 
part. Recently, a new Lindblad master equation adopted to the strong-coupling regime 
has been derived in Ref. |9] and it cures the unpredicted effects originated from the 
use of the ad hoc master equation (IHl) in the strong-coupling regime. The "improved^ 
Lindblad master equation can be written as [9]: 

p = -^[Hn,p]+J2^i'^l\3)(k\]p, (4) 

j,k>j 

where F^'^ = /t-^jy-^l (j| (a + o^) 1^) P are the relaxation coefficients [2U], \j) and cuj 
define the spectrum the Hamiltonian ([1]) i.e., \j) = Uj \j) and the eigenstates should 
be labeled according to increasing energy (label |j) such that Uk > coj for k > j). Also, 
note that as (7 — )■ 0, the standard dissipator can be obtained from Eq. (jlj). 

For the numerical solution of the master equation (|4]) with Rabi model, we 
have considered basis vectors of type lu)^^ = Ki^) where i,j = e,g and 
n = 0, 1, 2, . . . , d—1. The field Fock space is truncated for large enough d until the results 
converge [10]. In our calculations we have considered the resonant case, Uf = Ua = oj. 

Now, we investigate the dissipative dynamics of entanglement between one of the 
atom, (say A) and the cavity mode F in the strong-coupling regime. The reduced 
density matrix of A-F system can be obtained by taking a partial trace of p over the 
remaining part B: p^p = TtbP- In practice, the considered system has a dimension 
of (2 ® 00). Since we apply a convergence criteria to the field Fock space to solve the 
master equation (jl]), the atom-field system becomes a qubit-qudit (2 ® d) system. 

3. Atom-field entanglement in the ultrastrong-coupling regime 

Entangled states are vital for many quantum computation and communication protocols. 
Entanglement of formation (EOF) is one of the meaningful and physically motivated 
measure of entanglement. For two-qubit states, EOF is well known and has an 
elegant formula in terms of the entanglement monotone, concurrence [21]. For higher 
dimensional bipartite mixed states, except for some special cases 0, the evaluation 
of EOF is a formidable task which requires highly complex optimization process. 
Nonetheless, we can use recently proposed the lower bound of EOF [19] to estimate 

X One of the special case is the analytic calculation of EOF for a family of (2 ® d)-dimensional bipartite 
mixed states reduced from tripartite {2®2®d) pure states [17l [18] which will be discussed later in the 
text. 
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the amount of entanglement for the atom-field states in the strong-coupling regime. Its 
definition is based on the comparison of two strong separability criteria: The positive 
partial transpose and the realignment criteria. It was shown that EOF is tightly bounded 
for qubit-qudit system as [T9] : 





E{paf) > 



(A-i: 



where A = max 



[X 



if A = 1, 
if A G [1,2] 

1 - X)l0g2(l 



(5) 



X] 



and 



\\R{Paf)\\), H2{x) = -xlog2 
||G|| = Tr{GG'^)'2 is the trace norm. The matrix is the partial transpose with respect 
to the subsystem A. pj^ > is, in fact, sufficient for separability of bipartite (2 (g) 2) and 
(2 3) states. R{paf) is the realigned version of paf |l9] and is another operational 
criteria for separability which satisfies ||i?(pjj)|| < 1 for any bipartite separable state pij. 
These conditions state that if A > 1, the bipartite state is entangled. 




Figure 1. Time evolution of EOF for the qubits A and B (dashed line) and for the 
qubit A and the cavity mode F (solid line) for p(0) = l^eO) {geO\ initial state under 
RWA with g ~ OAuj and k = 0.2uj obtained by using the standard dissipator (|3]). 



We first investigate the dissipative dynamics of entanglement between the atom A 
and cavity mode F under RWA conditions. In Fig. [H we display the EOF between the 
atom A and cavity mode F as well as between the atoms A and B versus dimensionless 
time, ut for |\I/(0)) = \ge)^^ \0)p initial state with g = OAu and k = 0.2u by using the 
standard Lindblad master equation ([3]). One should note that for a single excitation 
initial state, the dimension of the whole Hilbert space will evolve to (2® 2® 2) due to the 
excitation number conservation of the RWA dynamics. Therefore, each party becomes 
a qubit-qubit system. For 2 qubit states, we can calculate EOF exactly [21] which 
is shown in Fig. [H Even from an initial factorized state, the atom-field interaction 
can create quantum correlations in time between each party. The peculiar behavior 
of correlations is seen in the long-time limit; the field becomes separable from the 
atoms {E{p^p) = E{p^^p) = 0), even the atomic state reaches a highly entangled 
steady state (E^p^^^) ^ 0.35). In fact, the behavior of the correlations in the steady 
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states has a simple explanation in terms of the decoherence-free subspaces (DFS) of 
the model; the atom-field system is driven by the cavity decay to a trapped mixed 
state involving only the two DFS states: the ground state, \ggO), and the dark state, 
|$~0) = -^{\egO) — \ geO)). The trapping properties are solely determined by the overlap 
between the atomic initial state and the subradiant state, |<I>~) (see Refs. [101 122] for 
details). More precisely, the population ($~| ^^^(0) |$~) is trapped in the atomic 
subradiant state, while the remaining part decays to the ground state of the atom-field 
system. Therefore, the atom-field steady states can be written as [lOJ [§|: 

p'' = [{1 - b) \gg) {gg\ + 6 |$-) ($- 1] ® |0) (0| , (6) 

where b = | p^s(O) |$^). Because of the presence of a maximally entangled 
subradiant state, |$~), the atoms can get entangled in the steady states. However, 
the cavity mode always becomes separable from the atoms in the steady states because 
of the separable nature of the field state in the two DFS. In fact, this description is 
only valid under RWA conditions and in the following we will show that the atom-field 
system can become entangled in the steady states under non-RWA conditions. 
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Figure 2. (a) Time evolution of lower bound of EOF between the qubit A and the 
cavity mode F for p{0) = \geO) {geO\ initial state in the strong-coupling regime with 
g = 0.25a; (sohd), g — OAuj (dashed), g = 0.6w (dotted) and k = 0.2w obtained 
by using the dissipator (|4]). (b) Ground state EOF between the qubit A and cavity 
mode F versus g/uj obtained from the monogamic relation (circles) and from the lower 
bound (rectangular). 



In Fig [21 we have investigated the dissipative dynamics of atom- field entanglement 
in the strong-coupling regime by using the dissipator, Eq. (jl]). Fig. [21^a) shows the 
dynamics of lower bound of EOF for the same initial state |^I/(0)) = \ge)^^ \0)p (Fig. [1]) 
and for several g/u values at k = 0.2a;. Similar to the RWA case (Fig. [T|), the atom-field 
system becomes entangled in a short time, but contrary to the RWA case, the field does 
not become separable from the atoms in the asymptotic time limit in the strong-coupling 
regime. In fact, in Ref. it was conjectured that for the extended Werner-like atomic 

§ Here we consider initial states that have no initial coherence between ground and dark states, i.e., 
{gg\ pab{0) \^^) — 0. Such states correspond to a wide range of initial states including X states [TO] , 
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initial states, the atom-field steady states have the same form of Eq. (jSj) with just l^'^'O) 



ggO), of the Rabi Hamiltonian: 



replaced by the ground state, 

p'' = {l-b) + 6 |$-0> ($-0| . (7) 

It is immediate to see that the atom-field entanglement, more precisely the role of CRTs, 
in the steady states stems from the ground state of Rabi Hamiltonian [23]. In Fig. Mjo), 
we have analyzed the entanglement property of the ground state for A-F system. One 
should note that the ground state forms a tripartite pure state of dimensions (2 ® 2 (gxi). 
It was shown that for a tripartite pure density matrix, EOF and quantum discord (or 
classical correlations [IT]) obey an important monogamic relation |18] : 

E{paf) = D{pab) + Sa\b, (8) 

where Sa\b = S{pab) — S{pb) is the von-Neumann entropy) and D{pab) is the 

quantum discord between the atoms which can be calculated analytically for the ground 
state [inilll]. Therefore, we can calculate ground state EOF between A and F exactly by 
using Eq. ([H]). The circle line in Fig. 2(b) displays the ground state EOF as a function of 
coupling strength g/uo for the resonance case obtained from monogamic relation, Eq. ([8]). 
EOF grows from zero (RWA regime) with the coupling up to ~ 0.5 and drops as g/u 
further increases after it approaches the maximum. In Fig. [2](b), we have also displayed 
the lower bound of EOF (rectangular line) between A and F for the ground state and 
it is found to be in agreement with the exact EOF [25] . 



4. Conclusions 



The non-RWA dynamics of atom-field entanglement has been investigated for a system 
composed of two qubits resonantly coupled to a single mode leaky cavity field. The 
strong atom-field interaction is found to induce atom-field entanglement in the steady 
states contrary to the weak-coupling one where RWA is valid. 
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